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We study temperature effect on anomalous viscosity of Graphene Hall fluid within quantum many-vortex 
hydrodynamics. The commonly observed filling fractions, u in the range 0 < < 2 is considered. An 

expression for anomalous viscosity dependent on a geometric parameter-Hall expansion coefficient, is obtained 
at finite temperatures. It arises from strained induced pseudo-magnetic field in addition to an anomalous 
term in vortex velocity, which is responsible for re-normalization of vortex-vortex interactions. We observed 
that both terms greatly modify the anomalous viscosity as well as an enhancement of weakly observed v 
fractions. Finite values of the expansion coefficient produce a constant and an infinite viscosities at varying 
temperatures. The infinities are identified as energy gaps and suggest temperatures at which new stable 
quantum hall filling fractions could be seen. This phenomenon is used to estimate energy gaps of already 
measured fractional quantum Hall states in Graphene. 
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I. INTRODUCTION 

Graphene is a monolayer of carbon atoms tightly 
packed into a two-dimensional (2D) honeycomb lattice. 
Since its discovery in 2004, graphene has attracted a 
great deal of attention mainly due to its exceptionally 
high crystal and electro nic q uality. Shear viscosity has 
been studied in graphen^^. The shear viscosity estab¬ 
lishes transverse velocity gradients that obstruct coherent 
electron flow. Under some conditions, one has relativis¬ 
tic particles in graphene forming a quantum Hall fluid; a 
strongly correlated sates of matter, which flows just like 
fluid and without such shearing resistance or dissipation. 
The viscosity measured in the quantum Hall fluid is what 
is termed Anomalous (also called Hall or dissipationless) 
viscositji^. Quantum fluids are particularly interesting 
especially due to the remarkable natural phenomenon 
occurring in superconductivity, superfluidity, ultracold 
atoms. Fractional Quantum Hall (FQH) Effect is yet 
another example. The ground states of FQH states are 
holomorphic in nature and gapped. These gapped states 
are characterized by a universal anomalous viscosity. The 
viscosity is revealed when stress tensor becomes sensitive 
to stress preserving deformations of the fluid. In this 
context, the origin of the anomalous term is due to fluid 
velocity diverging at microscopic scale which deforms a 
metric locally and causes dilatation of particle coordi¬ 
nates. The divergences due to individual particles are 
collectively manifested at macroscopic scale as an anoma¬ 
lous term. Another important term arising from metric 
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deformations at finite temperature is Hall expansion co¬ 
efficient. It is well established that graphene can respond 
to local deformations by producing strain which in turn 
induces giant pseudo-magnetic fields as much as 10 T in 
strained graphen^ and 300 T in graphene nanobubble^. 
The induced field allows vortices in the system to feel an 
effective magnetic field. It is the Hall expansion coeffi¬ 
cient that captures the contribution. 

Recently, there has been a great deal of interest and 
renewed focus on the anomalous viscosity, r]A of quan¬ 
tum fluids. In particular, re:P obtained a universal rela¬ 
tion TjA = {hp/4)S for FQH states of generic bulk sys¬ 
tems by studying the response to metric deformations. In 
Ref.l^r?^ (= (H/STr)^) was obtained for the FQH states 
of relativistic particles in graphene by electromagnetic 
and gravitational response. The same general result, in¬ 
cluding anomalous term, was also obtained by Abanov 
within effective hydro^namic theorjEl and within Euler 
vortex hydrodynamic^. However, all these results lack 
finite temperature effects which can have remarkable con¬ 
sequence on rjA and other transport properties. 

In the following, we consider flows of quantized Hall 
vortices in the graphene Hall fluid as elementary objects 
forming themselves a highly correlated quantum fluid. 
In the regime of long-wave slow motion, an accurate hy¬ 
drodynamic description becomes possible and does, in 
its own validity, not depend on the microscopic behav¬ 
ior of the electronic fluid. Vortex-vortex interactions will 
then be responsible for appearance of vortex EQH Ef¬ 
fect. The hydrodynamics of vortex matter presented 
here differs from Euler hydrodynamic^ by an anomalous 
term. Within this many-vortex approach, the anomalous 
term has been derived for bulk single component non¬ 
degenerate system^ without the Hall expansion coefR- 
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dent. 

The remaining of the paper is organized as follows. In 
Section [ll] we obtained vorticity velocity equation from 
Euler hydrodynamic equation through Boltzmann trans¬ 
port equation for graphene. We obtained a quantized 
Helmholtz-Kirchoff vortex solution. Vortex flux and mo¬ 
mentum conservation laws are subsequently derived. The 
stress tensor is deduced from which the anomalous vis¬ 
cosity is read off. In Section m we analyzed behavior 
of the anomalous viscosity under density and tempera¬ 
ture profiles for different filling fractions. We concluded 
in Section |III| highlighting possible applications of our 
results. 


of vortices is easily obtained from Eq. It is given by 

= + ( 4 ) 

i^j i ^ ' J ^ ' 

Equation Q is the Helmholtz-Kirchoff equation for vor¬ 
tices. We will specialize in the zero energy state of the 
system where K and K' components of the fluid decou¬ 
ple. Dynamics are then localized in either sublattice and 
the a dependence may be dropped. 

B. Temperature effect on dynamics of hall vortices and 
viscosity anomaly 


II. THEORETICAL MODEL 

A. Euler hydrodynamics and point vortices 

Two dimensional Euler hydrodynamics can be 
straightforwardly derived following Boltzmann transport 
equation at local equilibriurrfiSI, Particle distribution is 
reduced to continuity and Euler equations for density and 
momentum 

Dtp°‘ = 0 and Au“ -I- Vp“ = 0, (1) 

respectively. Where Dt = i9t -f u • V and u is macroscopic 
fluid velocity connected with the microscopic electron ve¬ 
locity, V = vp^/k. p is the partial pressure per density 
and a is the fluid component index. For graphene, a = 
(K't) Kf, K' I, K|). K,K' and are the valley and 
spin indexes, respectively. Taking the curl of the Euler 
equation Q, we get 


= 0. (2) 

Where the quantity w = V x u is vorticity which is non¬ 
zero for rotational incompressible flows, uj is considered 
as frozen into the fluid and constitute its own fluid. The 
continuity equation for the vortex fluid = 0, also 

holds. A solution for equation ([^ exists and consists of 
point-like vortices. The fluid velocity is 

u“ = -iOr+zy-(3) 

Where z = x + \y, u = Ux — iuy, T is the circulation 
and n which is identihed with cyclotron frequency of 
the fluid particles. Assuming a flow in which the cir¬ 
culation, r“ (= r“) is both minimal and chiral so that 
in the limit, W —>■ oo rotation is compensated by the 
large number of vortices. If magnetic field is present, the 
vortices of all components are smoothly distributed with 
fixed mean density, po = (I/7r)(D/r). Borhn-Summarfelt 
phase-space quantization leads to quantization of the cir¬ 
culation rrivE" = j3°‘h. Where is inertia vortex flow 
and I3°‘ is an integer. The equation for the slow motion 


In the quantum Hall regime of dissipationless flow, 
fluid particles do not carry heat flux, but vortices do. 
Vortices move in response to temperature gradient, VT. 
This can be captured in a vortex density dehned through 
momentum flux as P(r,T) = my{py{r,T),\}/2 


P(r,r) = ^{pv(r,T),v}. (5) 

Using the identity Tidlr) = d{l/z) and the ward identitjEl 
in Eq. we get 


P(r,r) = m^p^{r,T) 


u-f 


r 

4 


V* log p,{r,T) 


( 6 ) 


Equation ([^ is crucial in this studies. In particular, 
our results is based on the second term which dictates dis¬ 
cussions that follows. The term is responsible for anoma¬ 
lous behavior of the fluid when approaching a vortex. It 
is a quantum or micro-scale phenomenon which manifest 
itself at classical regime due to possible broken transla¬ 
tion symmetry associated with lattice scale deformations. 
Its presents renormalizes the vortex-vortex interactions 
in equation Q. It also creates stresses perpendicular to 
the fluid flow with no work nor dissipation. The asso¬ 
ciated transport coefficient (viscosity) is expected to be 
dissipationless. The momentum conservation following 
from broken translation invariance in the presence of ex¬ 
ternal forces yields 


Pa + V^Hafa — PvPa- (7) 


Where the stress tensor is H = 77ivPvU-u-|-77^(V*u-fVu*) 
V*(-) = V X (•). F = E + vxB is the external force. 
The kinetic coefficient rjA is the anomalous (dissipation¬ 
less) viscosity. Except for the temperature dependence, 
the dessipationless viscosity has the same structure as 
the universal relation obtained in literature for graphene 
Fractional quantum Hall stated, i.e, 

TD r 

rjA{r,T) = ^pAr,T). ( 8 ) 

In order to have a complete description of our system, 
one needs to quantize Eq. It leads to physical interpre¬ 
tation of the quantized circulation as the filling factor of 
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the fractional vortex Hall stated. A connection between 
r(= (fi/m^)S) and the so called shift, 5 of a Hall state 
is easily established. For Graphene, S = ~ l)/^- 

Except the re-definition of F, the quantized vortex veloc¬ 
ity assumes the same form as Eq. Taking the curl and 
and utilizing the relation V x u = 27rr(pv — Po), we have 


III. RESULTS, DISCUSSION AND CONCLUSION 

We now discuss the behavior of anomalous viscosity. 
Here, we observed how small temperature gradients 
affect viscosity of vortex fluid quantized on Hall states 
having filling fraction within Q < vq < 2. 


V X V = 27 rr 


p^{r, A*logpv(G T) 


(9) 


C. Effective magnetic field 


As we have pointed out in the introduction, geomet¬ 
ric deformations can lead to induction of giant synthetic 
magnetic fields. Vortices feel this contribution in addition 
to external magnetic field as an effective field, B^. The 
anomalous viscosity for the graphene-vortex system takes 
the shape, riA[r,t,T) = {eK/8Tr)By{r,t,T) with n = vqS 
and p = vcifiB/h). The Fourier transform of B^ in time- 
domain yields frequency-dependent viscosity 

VAir,uj,T) = ^B^{r,uj,T). ( 10 ) 

Where r = \zi — Zj\. To compute Bv{r,uj,T), we decom¬ 
pose it to B^{uj)B^{r,T)/B q. Bq is the applied field. 

First we need to compute B^ from Eqj^ To do that, 
we require that in the ground sate V x v = 0 which 
yields B^{r,T) = Bq- {H/4e){{l/vG)-2)A*\ogB^{r,T). 
Expanding B^ in small temperature gradients to get an 
iterative expression for the magnetic field, 


Hv(r,r) = Bo~ {AHogHv(r) + AHogJl - 7VT(r)] } 


( 11 ) 

Where ^ = 1/2 z/g —1 is the anomalous term. The Hall ex¬ 
pansion coefficient is defined as 7 = —{l/B^){dB^/dT). 
The minus sign is inferred from recent experimental re¬ 
sults on graphene at low-temperature^lSHAI. Because 
strain in graphene induces giant synthetic magnetic 
fields, at finite temperatures 7 can have great conse¬ 
quences on the system’s electronic transport properties. 
Thus, the coefficient can characterize stress and knowing 
it can be very critical in strain engineering. 

Finally, to obtain B^{ijj) we recall that the Magnus 
force acting on vortices has th form 


p^F = mvfli X (pqv - Pvu). (12) 


The term in brackets is written in such a way that the 
force stays constant in order that fluctuations in the zero 
energy states are bounded. Replacing v with equation 
^ and using p^F = P = —iiorriyp^v up to leading order 
in gradients expansion, we get 


Bo[l + {^/2){kef] 
(cu/0)2-[l + (e/2)(H)2]^' 


(13) 



Vc 



FIG. 1. Temperature gradient enhanced observations of weak 
fractions. To = 4.0 A, 7 = 0.9 u/fl = 0.1 and k£^ = 0.2. 
(top) 0 < I'G < 1 and (Bottom) 1 < vq < 2. 

In Figure 1 (a), the deviation of the 77 y 4 -curves from the 
mean viscosity, tja = measures the observability 

(strength) of a Hall state within 0 < i/q < 1- It is clear 
that both I'G = 1/3 and vq = 218 give the highest devia- 
tion s whi ch collaborate recent experimental results. See 
Ref^mm and references there in. Figure 1 (b) is similarly 
obtained for fractions within 1 < < 2 , except that the 

mean viscosity is rjA = —0.5577o.t- The negative viscosity 
may suggest why these fractions are weakly observed in 
experiments. 

In Figure 2 (a) and (b), we observe the anomalous 
viscosity over varying temperature, T. At some critical 
values, Tc, the viscosity grows to positive and negative 
infinities for vg > ^jl and vg > 4/7, respectively. Figure 
2 (c) and (d) was similarly obtained for fractions within 
1 < r-G < 2 . Above and below Tg, pa is robust (flat). 
This particular character places graphene as an impor¬ 
tant industrial material where temperature and viscosity 
are of utmost importance. 

In Figure ?? we plotted p in the neighborhood of in- 
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finites initially excluded from Figure We extracted the 
critical temperatures Tc giving rise to rjA infinities and 
study ??(T( 7 , vc) behavior as shown in Figure]^ The fol¬ 
lowing characterized the appearance of observed empty 
spaces in the plot. The wider empty spaces (contain¬ 
ing black squares) mean many fractions still remain hid¬ 
den which experiments have not yet measured. With im¬ 
proved precision and sensitivity, experiments should be 
able to observe these fractions. 

The intuitive meaning of the infinities in rjA can be 
clarified further. The infinities correspond to, up to 
a factor of some function, /(t'c) of the filling factor, 
the activation energy gap To estimate the energy 

gap, we write = f{vG)Tc and expand /(vg) in the 
neighborhood of one of the transition zones in Figure 
VG 1/2 and ~ 1. After applying some approxi¬ 
mations, we obtained /(vg) = — ^'g)- Where 

Ajyg is the degeneracy of z/G-fraction. For example, for 
i/G = 1/3, Ai /3 = 1 and A 1/3 = 18.14 AT. Similarly, 
one has A 2/5 = 5.47 AT, A 3/7 = 5.22 AT, A 4/7 = 3.95 AT, 
A 3/5 = 3.68 AT, A 2/3 = 9.25 AT. These values agreed fa¬ 
vorably with recent experimentJl^E^. 

In conclusion, we have computed dissipationless vis¬ 
cosity of quantum vortex Hall states in graphene within 
hydrodynamics using quantum many-vortex picture of 
Euler hydrodynamics. The hydrodynamics formalism al¬ 
lowed a great deal of simplifications as the microscopic 
theory is completely unnecessary and only few variables, 
p and V are employed. 

An important aspect of the dissipationless viscosity are 
the anomalous ^ and Hall 7 terms. Their effects come 
in two folds, (i) They reconstruct vortex density by in¬ 
troducing non-linearities in the fluid dynamics making 
flows at different FQHS distinguishable, (ii) ^ captures 
temperature variations due to strain effects through 7 
as soon as vortices are formed. The combined effect 
make viscosity robust over wider temperature range and 
can even enhanced observation of weaker FQH states. 
Controlling these parameters can enhance observability 
of weaker fraction and expose weaker ones. Knowing 
these parameters may be critical in strained-engineered 
devices. Since device sensitivity to stress-preserving de¬ 
formation depends on the factors. 


rjA contains infinities within some specific temperature 
ranges. These temperatures are associated with energy 
gaps. We estimated this energy gaps within 0 < < 2 

FQH states by analysing the temperature dependence of 
rjA- C and Hall 7 terms shift to higher values. We 
found Ai /3 = 18.14 AT that has been difficult to mea¬ 
sure or mostly underestimated in experiments. Our esti- 
mat ed gap s also compares favorable with literature. See 
] 7 efsiiSli£l and references there in. 

Finally, our results could be applied to strained- 
engineered devices to control viscosity. The studies can 
also guide future experiments towards observing new 
fractions. In particular, the temperature windows (white 
spaces in Figure]^ obtained may be probed, though away 
from the transition zones, for new fractions by controlling 
the Hall expansion coefficient parameter. Moreover, our 
work may resolve conflicts of different reported energy 
gaps. Specihcally, the A 1/3 gap. 
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FIG. 2. Normalized anomalous viscosity plots versus temperature at different fractions showing infinities corresponding to 
critical temperatures. 7 = 0.9 uj/Q. = 0.1 and kt = 0.2. (top) 0 < < 1 and (Buttom) 1 < j/g < 2. 
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FIG. 3. (Left) Normalized viscosity-temperature plot. Poles correspond to critical temperatures, Tc- (Right) j^A-Tc-r'G-plot. 
Full squares indicate critical temperatures. 7 = 0.9 cu/H = 0.1 and ki = 0.2. 

















































